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Abstract 

We calculate renormalization constants Zh,w, the Higgs and PT^-boson mass and tadpole 
counterterms in the on-mass-shell renormalization scheme to two loops in the heavy-Higgs- 
boson limit mn ^ Mw- Explicit analytic formulae are presented for the two-loop integrals 
with masses, which are not known in the literature. As an application, the analytic expression 
for the two-loop leading correction to the fcrmionic Higgs boson width is obtained. 



1. Introduction 



After the observation of the top quark signal at the Tevatron, the mechanism of the spontaneous 
electroweak symmetry breaking remains the last untested property of the Standard Model (SM). 
Although the recent global fits to the precision electroweak data from LEP, SLC and Tevatron seem 
to indicate a preference to a light Higgs boson uih = 149lg2^ GeV, uih < 450 GeV (95%C.L.) 
||^], it is definitely premature to exclude the heavy Higgs scenario. The reason is that a restrictive 
upper bound for niH is dominated by the result on A^ji, which differs significantly from the SAi 
predictions 0] . Without A upper bound on Mh becomes larger than 600 GeV , which is not 
far in the logarithmic scale from the value of the order of 1 TeV, where perturbation theory breaks 
down. In order to estimate the region of applicability of the perturbation theory, the leading two- 
loop 0{g^mff/M^^) electroweak corrections were under intense study recently. In particular, the 
high energy weak-boson scattering corrections to the heavy Higgs line shape Q, corrections to 
the partial widths of the Higgs boson decay to pairs of fermions [||, and intermediate vector bosons 
10, 1^ at two-loops have been calculated. All these calculations resort at least partly to numerical 
methods. Even for the two-loop renormalization constants in the Higgs scalar sector of the SAi 
IQ, complete analytical expressions are not known. In this paper we present our analytic results for 
these two-loop renormalization constants, evaluated in the on-mass-shell renormalization scheme 
in terms of transcendental functions ^(3) and the maximal value of the Clausen function Cl(7r/3). 



2. Lagrangian and renormalization 

The part of the SAi Lagrangian describing the Higgs scalar sector of the SAi in terms of the bare 
quantities is given by: 

- ^ [w>o + \4 + \hI + voHo + ^dv^y. (2.1) 

Here the tadpole counterterm 6v'^ is chosen in such a way, that a Higgs field vacuum expectation 
value is equal to vq 

»„ = ^ (2.2) 

g 

to all orders. 

Renormalized fields are given by 

Hq = ^fZ^ H, zo = ^/Z^ z. Wo = \/Z^w. (2.3) 

At two-loop approximation the wave function renormalization constants, tadpole and mass 
counterterms take the form 
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™2 _ 2 I 9^ , 2 (1) I r 2 (2) 

"iffo - + (167r2)2 ^ ■ 

Since weak coupling constant g is not renormalized at the leading order in the M^-boson mass 
counterterm is related to the Nambu-Goldstone wave function renormalization constant by the 
Ward identity Mwo = ZujMw- 

In the on-mass-shell renormalization scheme all the counterterms are fixed uniquely by the 
requirement that the pole position of the Higgs and I^-boson propagators coincide with their 
physical masses and the residue of the Higgs boson pole is normalized to unity. 

The one-loop counterterms equivalent to those used in 0, |^ are given by 




bm\^^^ _ m\i% { 3 Svr^/S / vr^ SvrVS 3^/3C 37r^/31og3 

.7(1) _ "^Ulf [3 , / 3 3vrV3 ^/3C 7r^/31og3 \] 



6Zm = -^\l 
w I 



Mw " 16 



Here the dimension of space-time is taken to be d = 4 + e and 

ei, = e-/^(ff)\ (2.6) 

In contrast to papers |9| we prefer not to keep the one-loop counterterms of 0{e) order, i. e. unlike 
in the conventional on-mass-shell scheme used in ^ , we require that the one-loop normalization 
conditions are fulfilled only in the limit e — > 0, where the counterterms of the order 0{e) do not 
contribute. Such a modified on-mass-shell scheme is equally consistent as the conventional one or 
the standard scheme of minimal dimensional renormalization, which assumes only the subtraction 
of pole terms at e = |]l^. (Moreover, in general one cannot subtract all the nonsingular 0{e) terms 
in the Laurent expansion in e, as they are not polynomial in external momenta.) The 0{e) one-loop 
counterterms considered in [^, |9| do can really combine with the 1/e terms at two-loop order to give 
finite contributions, but these contributions are completely canceled by the additional finite parts 
of the two-loop counterterms, fixed through the renormalization conditions in the on-mass-shell 
renormalization scheme. The reason is that after the inclusion of the one-loop counterterms all 
the sub divergences are canceled and only the overall divergence remains, which is to be canceled 
by the two-loop counterterms. The account of finite contributions coming from the combination 
of 0{e) one-loop counterterms with 1/e overall divergence just redefines the finite parts of the 
two-loop counterterms. An obvious advantage of this modified on-mass-shell scheme is that the 
lower loop counterterms once calculated could be directly used in higher loop calculations, while in 
the conventional on-mass-shell scheme for Z-loop calculation one needs to recalculate the one-loop 
counterterms to include all the terms up to C(e'~^), two-loop counterterms to include 0(e'~2) terms 
and so on. 
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3. Analytic integration 

The calculation of the Higgs and M^-boson (or Nambu-Goldstone w, z bosons) two-loop self en- 
ergies, needed to evaluate the renormalization constants ( |2.4| ), reduces to the evaluation of the 
two-loop massive scalar integrals 

J{k'^;niml,n2ml,n3ml,n4ml,n5ml) = J D^'^^PD^'^^Q (^P'^ -mj^ ' (3.1) 

xi^{P + kf-mlj i^{Q + kf-mlj [Q^-mlj i^{P - Qf - mlj 

and their derivatives J' at k'^ = m'j^ or at k'^ = 0. 

The most difficult is a calculation of the all-massive scalar master integral corresponding to the 
topology shown in the Fig. 1. 



Fig. 1. The two loop all-massive master graph. Solid line represents Higgs bosons. 



This integral has a discontinuity that is an elliptic integral, resulting from integration over the 
phase space of three massive particles, and is not expressible in terms of polylogarithms. However 
one can show |ll| that on-shell k"^ = m'jj or at the threshold k"^ = 9m^ this is not the case. We 
use the dispersive method ||l^, |l^ to evaluate this finite integral on the mass shell: 

m\ J{k^]m]j,m\,rn\,m]j,m]j) = aa{k^ /m\) + abik"^ /mjj), (3.2) 

where aa,b correspond to the dispersive integrals calculated, respectively, from the two- and three- 
particle discontinuities, which are itself reduced to one-dimensional integrals. The tanh~^ functions 



entering aa,b can be removed integrating by parts either in the dispersive integral |12|, or in the 
discontinuity integral |jl^ . By interchanging the order of integrations the latter representation gives 
the three-particle cut contribution ab as a single integral of logarithmic functions ||l^ . After some 
rather heavy analysis we obtain at fc^ 



^a(l) 




{y^ - 1) log(j/) _ vry 

y \/3 



^ C(3) - y vrC + vr^ log2 - - ^2 ^^^^^ 



dy 2 log 



y' + y + i 
y 



log(j/ + 1) ^ 3 y + 1) - 2 (y2 _ 1) log(y) ^^g^^) 



y 



y^ + y' + i 



y + i 



^C(3) + ^ ^ C - vrMog 2 + 3_ 



(3.3) 



(3.4) 
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Here C = Cl(7r/3) = Im lis (exp (ivr/S)) = 1.01494 16064 09653 62502 . . . 
As a result we find 

m]jJ{m]j;m\,m\,m]j,w?H,m\) = C(3)-^vrC 

= -0.92363 18265 19866 53695 . . . (3.5) 

The numerical value is in agreement with the one, calculated using the momentum expansion ||l4 
and with the numerical values given in 15 1. 



Given the value (|3.5D , the simplest way to calculate the derivative of the integral (|3.2| ) is to use 
Kotikov's method of differential equations |l^, ^ 

4t// 2 2 2 2 2 2\ ^ TT^ 

niH J {mjj;'mH,mH,mH,mH,mH) = 3^C'-C(3)-y 

= -0.17299 08847 12284 42069 . . . (3.6) 

All the other two-loop self energy scalar integrals contain "light" particles (Nambu-Goldstone 
or W, Z-bosons) and some of them are IR divergent in the limit Mw^z — > 0. In principle, one can 
calculate these IR divergent integrals in Landau gauge, where masses of Nambu-Goldstone bosons 
are equal to zero and IR divergences are represented as (double) poles at e = However, 
in order to have an additional check of the cancellation in the final answer of all the IR divergent 
log(M^2)-terms, we work in 't Hooft-Feynman gauge. For the infra-red finite integrals the correct 
answer in the leading order in rnjj^ is obtained just by setting Mw^z = 0. We agree with the results 
for these integrals given in |^ 0. The two-loop IR divergent integrals correspond to the topologies 
shown in the Fig. 2, which contain "massless" propagators squared. 



Fig. 2. The two-loop IR divergent graphs. Dashed line represents "light" particles. 



The relevant technique to handle these diagrams follows from the so-called asymptotic operation 
method [17|. According to the recipe of As-operation, the formal Taylor expansion in small mass 
M\y entering propagator (or its powers) should be accomplished by adding the terms, containing 
the J-function or its derivatives. The additional terms counterbalance the infra-red singularities, 
arising in the formal expansion of propagator. In our case we have 



1 1 M, 



+ 2—:^ + . 



2 

W 



[P^-M^f (P2)2 (p2)3 

+ Ci{Mw)5^''\P) + C2{Mw)dH^''\P) + ... (3.7) 
Here the first coefficient functions Ci{Mw) read 
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This equality is to be understood in the fohowing sense |17]. One should integrate both parts of the 
equation multiplied by a test function and then take the limit d — > 4. If one keeps in the expansion 
all terms up to order M^, the resulting expression will represent a correct expansion of the initial 
integral to order o(M^). To obtain the leading contribution to the diagrams Fig. 2, it suffices 
just to take the first term of the Taylor expansion and, correspondingly, the first "counterterm" 
Ci5(^)(P). 

Finally, the combination of Mellin-Barnes representation and Kotikov's method gives the fol- 
lowing answer, corresponding to the first graph in Fig. 2 neglecting the terms of order 0(M^/m^): 

m^J(m^;2M^,M^,0,M^,m^) = \ 1^ - '-^ + ^ log{^) (3.9) 

1 




m 



H 



m\j J' {ml- 2 Ml,, Mir, Ml,, ml) = ( + 2z tt - ^ + log(^) ) (3.10) 

Ml, 1 ^^Ml,. 
+ log^!-) - - log ' 
mjj 

The integral diverges as 1/e, while if we would set Mw = from the very beginning, it would 
diverge as The integral corresponding to the second graph in the Fig. 2 up to 0{Ml, /ml) is 

JKh;2A4,,0,0,a4,,„4) = {g^l + i('2 1og(:^)-l") (3.11) 

2 12 ^^ml^ 2 ^ ^ml 'j 

The two-loop vacuum integrals needed to evaluate the tadpole counterterm 6v'^^'^^ have been 
calculated in IT 




4. Results 

The analytic results for the two-loop renormalization constants are: 

= TTHTf — + 162 - 3 7r2 + 60y3C7 ; (4.1) 

Id Myj/ y e e j 

5ml^'^^ _ mf^^|//576 144 vr ^/3 - 1014 
ml ~ 64 1^12" + I 

99 

+y - 252C(3) + 877r2 -2197rV3 

+156 7rC + 204\/3C^; (4.2) 

6ZP = ^f^-I^-126C(3) 
^ 64M^ ^ 
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(2) 
W 



Mw 



25 vr' 

5Z^) - 



76 7r\/3 + 787rC + 108V3C ; 



64M|; U " 8 " Y + - 6 V3C I . 



(4.3) 
(4.4) 



For comparison we have also calculated these counterterms following the renormalization scheme 
iQ, |5| and keeping the 0{e) terms and found complete agreement with their (partly numerical) 



results. In this scheme the renormalization constants (4.1), ( f4.2[ ) look a bit more complicated due 
to the presence of the additional 7r\/31og3 terms 



5v 



2(2) 



5m 



2 (2) 
H 



m 



H 



m%CH ( 367r ^3-84 



16 M, 



+ 90 - 12 vr^ - 12 v^C 



w 



-367r\/3 + 18 7r\/3 log3 ; 



m 



jj^H [ 576 1447r\/3- 1014 



+ 



+- 



2439 



252C(3) + 637r2 - 363 7ry3 



(4.5) 



+156 7rC-84V3C + 72 7rV3 log3 . 



(4.6) 



The wave function renormalization constants Zh,w are identical in these two schemes. 

As an example of the physical quantity, for which all the schemes should give the same result, 
we consider the two-loop heavy Higgs correction to the fermionic Higgs width ||5|, ^ . The correction 
is given by the ratio 



Zh 



1 + 2 



+- 
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5M, 



{!)■ 
W 



Mw 



2\2 



16 vr 



+2 5Z^^^ 



5M, 



(1) 
w 



6M, 



(1) 
w 



Mw \ Mw 



5Z 



(1) 

H 



5M, 



(1)- 
w 



Mw 



6M, 



(2) 
W 



Mw 



Substituting (U), (gj 
Zh 



1 + 



we find 

1 g' 



m. 



16 V16 7r2M2 



+ 



w 



^ I- ( 13 - 27r\/3 
8167r2M2,V 

^3 - 63 C(3) - + ^ + 39 . C + 57 ^/3 . 



(4.7) 



Again, we find complete agreement with the numerical result H] and exact agreement with the 
result 0, taking into account that their numeric constant is just minus our integral (|37 
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